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Type II superconductors, consisting of superconducting domains embedded in a normal or in- 
sulating matrix, undergo a rounded phase transition. Indeed, the correlation length cannot grow 
beyond the spatial extent of the domains. Accordingly, the thermodynamic properties will exhibit 
a finite size effect. It is shown that the specific heat and penetration depth data of a variety of type 
II superconductors, including cuprates, exhibit the characteristic properties of a finite size effect, 
arising from domains with nanoscale extent. The finite size scaling analysis reveals essential features 
of the mechanism. Transition temperature and superfluidity increase with reduced domain size. The 
combined finite size and isotope effects uncover the relevance of local lattice distortions. 
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I. INTRODUCTION 

Since the discovery of superconductivity in cuprates by Bednorz and Miiller [1] a tremendous amount of work has 
been devoted to their characterization. The issue of inhomogeneities and their characterization is essential for several 
reasons, including: First, if inhomogeneity is an intrinsic property, a re-interpretation of experiments, measuring an 
average of the electronic properties, is unavoidable. Second, inhomogeneity may point to a microscopic phase separa- 
tion, i.e. superconducting grains, embedded in a non-superconducting matrix. Third, there is neutron spectroscopic 
evidence for nanoscale cluster formation and percolative superconductivity in various cuprates [2,3]. Fourth, nanoscale 
spatial variations in the electronic characteristics have been observed in underdoped Bi 2 Sr 2 CaCu 2 08 + a with scanning 
tunnelling microscopy (STM) [4-7]. They reveal a spatial segregation of the electronic structure into 3nm diameter 
superconducting domains in an electronically distinct background. On the contrary, a large degree of homogeneity 
has been observed by Renner and Fischer [8]. As STM is a surface probe the relevance of these observations for 
bulk and thermodynamic properties remains to be clarified. Fifth, in YBa2Cu307_,5, MgB2, 2H-NbSe2 and Nb77Zr23 
considerably larger grains have been established. The magnetic field induced finite size effect revealed lower bounds 
ranging from L = 182 to 814A [9]. Sixth, since the change of the lattice parameters upon oxygen isotope exchange 
is negligibly small, the occurrence of a significant change in the inhomogeneities spatial extent, will provide clear 
evidence for the relevance of local lattice distortions in mediating superconductivity [10]. 

Here we review the attempts to probe the granularity of type II superconductors by means of thermal fluctuations 
and magnetic fields [9,10]. It is well-known that systems of finite extent, i.e. isolated superconducting grains, undergo 
a rounded and smooth phase transition [11]. As in an infinite and homogeneous system the transition temperature 
T c is approached the correlation length £ increases strongly and diverges at T c . However, when superconductivity 
is restricted to grains with length scale Lj in direction i, £j cannot grow beyond L.;. In type II superconductors, 
exposed to a magnetic field Hi, there is an additional limiting length scale = y$o/ (fiHi) with a 3.12, related 
to the average distance between vortex lines [9]. Indeed, as the magnetic field increases, the density of vortex lines 
becomes greater, but this cannot continue indefinitely, the limit is roughly set on the proximity of vortex lines by the 
overlapping of their cores. Due to these limiting length scales the phase transition is rounded and occurs smoothly. 
Indeed, approaching T c from below the transverse correlation length Q in direction i and its real space counterpart 

£T = y^CjCfei where i ^ j ^ k, increase strongly. However, due to the limiting length scales Li and L^ i — ^/$o/ (a Hi), 
it is bounded and cannot grow beyond 

er (*p) = & \t P i ~^Li, 

y/ST (*p) = \f&^ 1 \hr = y/*o/ (aH k ) = L Hk , i ? j + k, 

where t p = 1 — T p /T c and v denotes the critical exponent of the correlation lengths. Beyond the mean-field approx- 
imation it differs from v — 1/2 and a w 3.12 is a universal constant [9]. In this context it is important to recognize 
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that the confinement effect of the magnetic field in direction i on fluctuations within a region Ln i acts only in the 
plane perpendicular to the field. As a remnant of the singularity at T c the thermodynamic quantities exhibit a so 
called finite size effect, i.e. a maximum or an inflection point at T p . Two limiting regimes, characterized by 



L «. = \ h± = KW-'^Vl , (2) 



can be distinguished. For Ln i < LjLk the magnetic field induced finite size effect limits y ^ £ fc to grow beyond 

L,Hi, while for > yjLjL^ the superconducting grains set the limiting length scales. Since Ln i can be tuned 
by the strength of the applied magnetic field, both limits are experimentally accessible. Ln i < y/LjLk is satisfied 
for sufficiently high and > yJhjL^ for low magnetic fields. Thus, the occurrence of a magnetic field induced 
finite size effect requires that the magnetic field and the length scales of the superconducting grains satisfy the lower 
bound Hi > $ / (aLjLk). Since superconductors fall in the experimentally accessible critical regime into the 3D-XY 
universality class [25], with known critical exponents and critical amplitude combinations, we take these properties 
for granted [12]. They include the exponents 

a = 2 - iv = -0.013, v = 0.671 (3) 

and the critical amplitude combinations 



where 



W k > A±v ^ = ^ > ^ = (4) 

A + 

— = 1.07, R~ = 0.815 (5) 

a and are the critical amplitude of the specific heat singularity, c = (A^/a) \t\~ a + i? ± , while £^ and Aoi are the 
critical amplitude of the correlation length and penetration depth in direction i. In the homogeneous system these 
length scales diverge as £7" = ^ \t\~", \ = Xoi \t\~ u ^ 2 and V~ is the correlation volume. 



II. SPECIFIC HEAT 



In Fig. la we displayed the data of Roulin et al. [13] for the temperature dependence of the heat coefficient of a 
YBa 2 Cu 3 06.9 single crystal with T c = 92.56K at various magnetic fields applied parallel to the c-axis (H c ). As a 
remnant of the zero field singularity, there is for fixed field strength a broad peak adopting its maximum at T p which 
is located below T c . As T p approaches T c , the peak becomes sharper with decreasing H and evolves smoothly to the 
zero-field singularity, smeared by the inhomogeneity induced finite size effect. Since T p decreases systematically with 
increasing field, condition (2) applies in the form Lh c < V 'L a Lt,and there is a magnetic field induced finite size effect. 
As the data of Roulin et al. [13] are rather dense and extend close to zero field criticality, a detailed finite size scaling 
analysis appears to be feasible. In Fig. lb we plotted the data in terms of (c (T, H c ) /T — B~) / (^A~ \t\~ a /aj versus 

t/H^ 2v . Apparently the data falls on a single curve, which is the finite size scaling function g (y) defined in terms of 

- B") / (A- \t\- a ) =g(y),y = t ($„/ (aH c fej 2 ) ) ^ . (6) 

A finite size effect requires that g(— 1) = 1 at T p , where y = y p ~ — 1, while g (±y — > 0) oc |j/P' a ' holds for a < and 
g (±y — > oo ) = 1. To qualify the remarkable collapse of the data we note that the estimates for A~ = A~ /a and B~ 
have been derived from the zero field data displayed in Fig. la. The solid line for t < is c/T = A~ \t\~ a + _B~with 
A- = A-/T c a = -0.073(J/K 2 gat), B~ = 0.179(J/K 2 gat), B+ = 0.181(J/K 2 gat), T c = 92.56K and the critical 
exponent a = —0.013, to indicate the inhomogeneity induced deviations from the leading zero field critical behavior of 
perfect YBa2Cu306.g. Since the data collapse onto the finite scaling function has been achieved without any arbitrary 
fitting parameter, the existence of the magnetic field induced finite size effect is well confirmed. Furthermore, Fig. la 
shows how the broad anomaly in the specific heat coefficient sharpens, while the maximum height at T p (H c ) increases 
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with reduced field strength, evolving smoothly to the zero field peak, rounded by inhomogeneities. From the finite 



size scaling plot displayed in Fig. lb we also deduce that with £ Qah = 14. 4A, that at y p , where |t| /if, 
scaling function is close to 1, as required. Furthermore, the solid line which is (c(T,H c ) /T — B 

0.957 (— t/H^ 2 \ , confirms the divergence of the scaling function in the limit ±y — > 0. 
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FIG. 1. (a): Specific heat coefficient c/T versus t of YBa2Cu30e.9 with T c = 92.56 K at various magnetic fields derived 
from the data of Roulin et al. [13]. The solid lines are c/T = A \t\~" + B with a and A + /A~ given by Eqs.(3) and the 
parameters listed in the text. The deviations of the zero field data from these lines reveal around t = the finite size effect 

,l/2w 



due to inhomogeneities; (b): ( c (T, H c ) /T - B - ) / (k~ \t\ a ) versus t/H?" for >T>;i_>< 'u :1 ( ) (i ., wit!) T, 



from the data of Roulin et al. [13] with H in T. The solid line is f c (T, if c ) /T - B~ 
a = —0.013, confirming the divergence of the scaling function in the limit ±y — > 0. 
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The shift of T p with increasing H c is according to Eq.(l) given by 

l/2v 



t m = 



(7) 
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FIG. 2. |t p | versus H c (•) and H ab (a) for YBa 2 Cu 3 06.9 with T c = 92.56K derived from the data of Roulin et al. [13]. The 
solid lines are \t p {H ab )\ = 0.0033 if^ 2 " and \t p {H c )\ = 0.0132 H] /2v with v listed in Eq.(3). 



In Fig. 2 we show t p versus if a h and ff c for YBa2Cu30g.9 single crystal with T c — 92.56K derived from the data of 
Roulin et al. [13]. The solid lines are \t p (H ab )\ = 0.0033 H^ 2 " and \t p (if c )| = 0.0132 Hl /2v with if, revealing that 
the critical regime is attained. With v given by Eq.(3) and a = 3.12 these fits yield for the critical amplitudes the 
estimates 
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«o~c = 5.7A, Co" b = 14.4A, V- = 2 = 469 A 3 . (8) 

To estimate the correlation volume from the specific heat coefficient in terms of the universal relation, A V~ = 
(i? ± ) 3 (Eq.(4) we note that A~ = T c a A~ , A-{cm- 3 )= (lO 7 /k B /V gat ) A" (mJ/K/cm 3 ), V gat = 8cm 3 and A~ = 
-0.073(J/K 2 gat), corresponding to the solid line in Fig. fa, give A~ = 7.96 10~ 4 A -3 , yielding with R~ = 0.815 
(Eq.(5)) the correlation volume V~ « 680 A 3 , which is reasonably close to the value derived from the magnetic field 
induced finite size effect (Eq.(8)). Since at the lowest attained fields H c — 0.25T and H a h = 2T a shift from T c to 
T p is still present (see Fig. la), we derive from Eq.(2) for the length scales of the sample inhomogeneities the lower 
bounds 



L ab > 515A, ^L ab L c > i82A. (9) 

As in zero field the correlation volume V~ |tp| — " cannot grow beyond the volume Vi of the inhomogeneities, this 
limiting volume scale is obtained from Vi = V~ \t p \ 3v and t p evaluated in zero field. With V~ = 469A 3 and 
\tp\ = 0.0025 taken from Fig. la, we obtain the estimates Vi — 7.5 10 7 A 3 and V^ 3 = 422A which is consistent with 
the aforementioned lower bounds. 

Next we turn to Nb77Zr23 [14], 2H-NbSe2 [15] and MgB2 [16], type II superconductors supposed to have comparative 
large correlation volumes. In Fig. 3 we displayed \t p \ versus H derived from the respective experimental data. In 
contrast to the corresponding plot for nearly optimally YBa2Cu30y_5 (Fig. 2), the data points to a linear relationship. 
Consequently, the critical regime, where \t p \ oc H l / 2v with v w 2/3 holds (Eq.(7)), is not attained. Nevertheless, 
there is clear evidence for a magnetic field induced finite size effect, because T p shifts monotonically to a lower value 
with increasing magnetic field. Since the data points to an effective critical exponent v m 1/2 which applies over 
an unexpectedly extended range, we use Eq.(7) with v = 1/2 to derive estimates for the amplitude of the respective 
transverse correlation lengths and the correlation volume in terms of 



with a = 3.12. 




H(T) 

FIG. 3. \t p \ versus H for Nt>77Zr23 [14], 2H-NbSe2 [15] and MgB2 [16] derived from the respective references. The straight 
lines are Eq.(10) with the parameters listed in Table I. 

The respective straight lines in Fig. 3 are this relation with the parameters listed in Table I. For comparison we 
included the corresponding parameters for YBa2Cu30g. 9 [13], YBa2Cu306.6i (K,Ba)Bi03 [17] and HgBa2Ca2Cu307_5 
[19]where the critical regime is attained. It is evident that Nb77Zr23, 2H-NbSe2 and MgB2 are type II superconductors 
with comparatively large correlation lengths. Compared to YBa2Cri306.9 and YBa2Cu306.6 the correlation volume 
is 3 orders of magnitude larger, which renders the amplitude of the specific heat singularity very weak (see Eq.(4). 
Nevertheless, the unambiguous evidence for the magnetic field induced finite size effect reveals that fluctuations, even 
though not critical, are at work. For this reason there is no critical line T C 2 (H) of continuous phase transitions, 
but a line T p (H) where the specific heat peak, broadening and decreasing with increasing field, adopts its maximum 
value and the correlation length attains the limiting magnetic length scale Lh^ Because the fluctuations are also 
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subject to the finite size effect arising from inhomogeneities with length scale Lj, the magnetic finite size effect is 
observable as long as yjLjLk > Ljj ( ■ The resulting lower bounds for the length scale of inhomogeneities affecting the 
thermodynamic properties are also included in Table I. Noting that these bounds stem from studies where no attempt 
was made to explore the low field behavior, it is conceivable that the actual length scale of the inhomogeneities is 
much larger. To our best knowledge, the only absolute reference stems from the finite size scaling analysis of the 
zero field specific heat data of nearly optimally doped YBa 2 Cu307_,5, where L was found to range from 290 to 419A 
[25,21]. Interestingly enough, the largest bound found here applies to the cubic superconducting alloy Nb77Zr23. 
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< 69 



Table I: Summary of the estimates for the critical amplitudes of the correlation lengths, correlation volume V c — 
i^Oab) £oc' anisotropy 7 = £o~afc/£cic' lower bounds for the length scales L ab and \JL ab L c of inhomogeneities, derived 
from experimental data by means of the magnetic field induced finite size effect. For YBa 2 Cu30g.6 we used 7 = 20, 
taken from Janossy et al. [20] to estimate ^Q ab . 

On the contrary, the specific heat measurements reveal that in the Bismuth- and Thallium based cuprates T p appears 
not to shift up to 14T, for fields applied parallel or perpendicular to the c-axis [18,22,23]. To derive detailed estimates 
we consider the single crystal data for the specific heat coefficient of Junod et al. [23] for Bi 2 .i 2 Sri.7iCai. 22 Cui.g508+5- 
To reduce artifacts of the huge background we displayed in Fig. 4 the data in terms of (c (H, T) — c (H = 0, <)) /T versus 
T for fields applied parallel to the c-axis. Apparently, there is no shift of the peak around T p — 85. IK with increasing 
magnetic field up to H — 14T. Noting that the same behavior was found for fields perpendicular to the c-axis [23], 
these data lead with inequality (2) to the upper bounds 

L ab < 69A, ^L ab L c < 69A, (11) 

uncovering nanoscale superconducting grains, consistent with the length scale of the inhomogeneities observed with 
STM spectroscopy [4-7]. As STM is a surface probe, our analysis establishes that these grains are not merely an 
artefact of the surface, but a bulk property with spatial extent, giving rise to a rounded thermodynamic phase 
transition which occurs smoothly. 



0.0002 



^ -0.0002 
t; -0.0004 

P -0.0006 

o" 

I -0.0008 

X -0.0010 
I— 

I -0.0012 

o 

-0.0014 




76 78 80 82 84 86 
T(K) 



90 92 94 



FIG. 4. (c(H c ,T) — c(H c — 0,T)) /T versus T of a Bi2.i2Sri.7iCai.22Cui.950g+5 single crystal taken from Junod et al. 
[23]. The magnetic field is parallel to the c-axis 
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III. PENETRATION DEPTH 



Considering again the 3D-XY critical point, extended to the anisotropic case, the penetration depths and transverse 
correlation lengths in directions i and j are universally related by [25,26] 

I = 16 ^ T . (12) 

When the superconductor is inhomogeneous, consisting of superconducting grains with length scales Li, embedded 
in a non-superconducting medium, £* does not diverge but is bounded by = (£7/) < L\, where i ^ j ' ^ k. The 
resulting finite size effect is clearly seen in the microwave surface impedance data for X 2 b (T = 0) /X 2 b (T) versus T 
of Jacobs et al. [24], displayed in Fig. 5a. The solid curve indicates the leading critical behavior of the homogeneous 
system. A characteristic feature of this finite size effect is the occurrence of an inflection point at T p w 87K, giving 
rise to an cxtrcmum in d (\ 2 ab (T = 0) /X 2 ab (T)) jdT. Here Eq.(12) reduces to 

1 1 _ Wir 3 k B T p 

X 2 ab (T p )~ X a (T p )X b (T p ) <!> 2 L C ■ 1 j 

With X ab (T = 0) = 1800A as obtained from ^SR measurements [28], T p w 87K and X 2 b (T = 0) /X 2 ab (T P ) = 0.066 
we find L c w 68A, consistent with lower bounds (11) derived from the specific heat data. Thus, although the 
superconducting grains are of nanoscale only, a very small correlation volume makes the critical regime attainable. 
Clear evidence for 3D-XY critical behavior was also observed in epitaxially grown Bi2212 films [27]. 




FIG. 5. (a)Microwave surface impedance data for X 2 ab (0) /\ 2 ab (T) (O) and d (\ 2 ab (0) /X 2 ab (T)) jdT (□) versus T of a 
high-quality Bi 2 Sr 2 CaCu208 + i single crystal taken from Jacobs et al. [24]. The solid line is \ 2 ab (0) /X 2 ab (T) = 1.2 (1 - T/T c ) 2/3 
and the dash-dot line its derivative with T c = 87. 5K, indicating the leading critical behavior of the homogeneous system. The 
dotted line is the tangent to the inflection point at T p fa 87K, where d (\ 2 ab (T = 0) /X 2 b (T)) /dT is maximum; (b) Finite size 
scaling function g(y) — (\o a b / X ab (T)) 2 \t\^ 1 ' versus y = t/\t p \ for the data shown in Fig. 5a. The solid line in the inset is 
Eq.(15) with g = 1.6. 



To substantiate the occurrence of a finite size effect we explore the consistency with the associated scaling function. 
In the present case it is defined in terms of 

( a^i%) 2 '' r = 9 {y) ' y = szgn {t) ltl (I 1 ) UV = sign {t) f p (14) 

For t small and L c — > 00, so that ±y — > 00 it should tend to g (y — > 00) = land <? (j/ — > —00) = 0, respectively, while 
for £ = and L 7^ it diverges as 

5(y-o)= ffo y- I/ = flo(^) , (15) 
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so that in this limit (Xoab/^ab (T C ,L)) = go \tp\ V = .9o?oc/^c- As expected, a sharp superconductor to normal state 
transition requires domains of infinite extent. Moreover at t p , y = 1 and d (Xoab/^ab (T, L)) 2 jdt = 0. Accordingly, 
there is an inflection point at t p . Since the scaling function g (y) depends on the type of confining geometry and on 
the conditions imposed (or not, in the case of free boundaries) at the boundaries of the domains, this applies to the 
amplitude go as well. In Fig. 5b we displayed the finite size scaling resulting from the data shown in Fig. 5a. In view of 
the fact that the data satisfies the limiting behavior of the finite size scaling function remarkably well, there is strong 
evidence for a finite size effect. 

IV. FINITE SIZE AND ISOTOPE EFFECTS 

Recently we explored the effect of oxygen isotope exchange in Yi_ x Pr x Ba2Cu307_,5 on the inhomogeneity induced 
finite size effect by means of in-plane penetration depth measurements [10]. Defining the relative oxygen isotope shift 
of a physical quantity X as AX/X = ( 18 X — 16 X)/ 16 X) we note that the shifts are not independent but according 
to Eq.(13) related by 

AL C = AT Pc | AA^T^ (ig) 
L c T Pa X ab (T Pc ) 

From the resulting estimates, listed in Table II The resulting estimates are summarized several observations emerge. 
First, L c increases systematically with reduced T Pc . Second, L c grows with increasing x and upon isotope exchange 
( 16 0, 18 0). Third, the relative shift of T Pc is very small. This reflects the fact that the change of L c is essentially 
due to the superfluid, probed in terms of \ 2 ab . Accordingly, AL C /L C w AX 2 b /X 2 b for x = 0, 0.2 and 0.3. Indeed the 
relative shifts of T Pc , X 2 b (T Pc ) and L c are not independent. 



X 





0.2 


0.3 


&T Pc /T Pc 


-0.000(2) 


-0.015(3) 


-0.021(5) 


AL Pa /L Pc 
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ib T Pc (K) 
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66.0(2) 
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16 V (A) 
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1820(20) 


2310(30) 



Table II: Finite size estimates for le T Pc , ls T Pc , 16 X 2 ab ( 16 T Pc ) / le X 2 ab (0) and 18 X 2 ab ( ls T Pc ) / le X 2 ab (0), and the result- 
ing relative shifts AT Pc /T p &nA AX 2 ab (T Pc ) /X 2 b (T Pc ) for an content of 89%. XQ L p J xii L Pc and AL Pc /L Pc follow 
from Eq.(13). 16 A a6 (0) are ^SR estimates [29] 

To appreciate the implications of these estimates, we note that for fixed Pr concentration the lattice parameters 
remain essentially unaffected [30,31]. Accordingly, an electronic mechanism, without coupling to local lattice distor- 
tions, implies AL C — 0. On the contrary a significant change of L Pc upon oxygen exchange uncovers the coupling to 
local lattice distortions involving the oxygen lattice degrees of freedom. A glance to Table I shows that the relative 
change of the domains along the c-axis upon oxygen isotope exchange is significant, ranging from 12 to 16%, while the 
relative change of the inflection point or the transition temperature is an order of magnitude smaller. For this reason 
the significant relative change of L c at fixed Pr concentration is accompanied by essentially the same relative change 
of X 2 b , which probes the superfluid. This uncovers unambiguously the existence and relevance of the coupling between 
the superfluid, lattice distortions involving the oxygen lattice degrees of freedom. Potential candidates are the Cu-0 
bond-stretching-type phonons showing temperature dependence, which parallels that of the superconductive order 
parameter [32] . Independent evidence for the shrinkage of limiting length scales upon isotope exchange stems from 
the behavior close to the quantum superconductor to insulator transition where T c vanishes [33]. Here the cuprates 
become essentially two dimensional and correspond to a stack of independent slabs of thickness d s [34,35]. It was 
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found that the relative shift Ad s /d s upon isotope exchange adopts a rather unique value, namely Ad s /d s w 0.03 [33]. 
Although the majority opinion on the mechanism of superconductivity in the cuprates is that it occurs via a purely 
electronic mechanism involving spin excitations, and lattice degrees of freedom are supposed to be irrelevant, the rela- 
tive isotope shift AL C / 16 L C w AA„ fc / 16 A„ fc w 0.15 uncovers clearly the existence and relevance of the coupling between 
the superfluid, lattice distortions and anharmonic phonons which involve the oxygen lattice degrees of freedom. 

In contrast to YBa2Cu 3 07_5, (K,Ba)BiC"3, MgB 2 , 2H-NbSe2 and Nb77Zr 2 3, where the lower bounds for the length 
scale L of the superconducting grains ranges from 182 A to 814 A, we have seen that the data for Bi2Sr2CaCu20s+5 
single crystals and Yi_ x Pr x Ba2Cu307_,5 uncovers the existence of nanoscale inhomogeneities, which are not merely 
an artefact of the surface, but a bulk property with spatial extent, giving rise to finite size effects and with that to 
a rounded thermodynamic phase transition which occurs smoothly. While there are many questions to be answered 
about the intrinsic or extrinsic origin of the inhomogeneities, the existence and the nature of a macroscopic super- 
conducting state, we established that a finite size scaling analysis yields the spatial extent of the inhomogeneities. 
Furthermore, we established that the change of the spatial extent of the inhomogeneities upon oxygen isotope exchange 
uncovers the relevance of lattice degrees of freedom in mediating superconductivity. 
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